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ABSTRACT 

Transitional circumstellar disks around young stellar objects have a distinc- 
tive infrared deficit around 10 microns in their Spectral Energy Distributions 
(SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggest- 
ing dust depletion in the inner regions. These disks have been confirmed to have 
giant central cavities by imaging of the submillimeter (sub-mm) continuum emis- 
sion using the Submillimeter Array (SMA). liowever, the polarized near-infrared 
scattered light images for most objects in a systematic IRS/SMA cross sample, 
obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, 
in clear contrast with SMA and Spitzer observations. Radiative transfer mod- 
eling indicates that many of these scattered light images are consistent with a 
smooth spatial distribution for micron-sized grains, with little discontinuity in 
the surface density of the micron-sized grains at the cavity edge. Here we present 
a generic disk model that can simultaneously account for the general features 
in IRS, SMA, and Subaru observations. Particularly, the scattered light images 
for this model are computed, which agree with the general trend seen in Subaru 
data. Decoupling between the spatial distributions of the micron-sized dust and 
mm-sized dust inside the cavity is suggested by the model, which, if confirmed, 
necessitates a mechanism, such as dust filtration, for differentiating the small 
and big dust in the cavity clearing process. Our model also suggests an inwardly 
increasing gas-to-dust-ratio in the inner disk, and different spatial distributions 
for the small dust inside and outside the cavity, echoing the predictions in grain 
coagulation and growth models. 

Subject headings: protoplanetary disks — stars: pre-main sequence — circum- 
stellar matter 
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Introduction 



Transitional disks have a small or no excess from ~l/im to ^lO/xm relative to thei r 
full disk cousins, but a significant excess at longer wavelength (jWilliams fc Ciezal 1201 ll ). 
suggesting cleared out i nner disk. This iii t erpretation dates bac k to the era of the Infrared 
Astronomical Satellite (IStrom et al.l Il989l ISkrutskie et al.lll990l ). and was later developed 
with the help of detailed near-infrared (NIR) to mid-infrared (MIR) spectra provided by the 
Infrared Spectrograph (IRS) on-board Spitzer Space Telescope. Detailed radiative transfer 
modeling suggests that this kind of SED is consistent with disk models which harbor a central 
(partially) depleted region (i.e. a cavity or a gap), while a "wall-like" structure at the outer 



edge of this region can be responsible for the abrupt rise of the SED at MIR ( iCalvet et al. 
2OO5I : lEspaillat et al.l 120071 1^. 



The disk -(-cavity model based on SED-only fitting usually contains large uncertainties, 
because the SED samples the emission from the whole disk; by tuning the ingredients in 
the fitting, one could fit th e I RS SED with differen t models (see the example of UX Tau A, 
Espaillat et al.ll2010l . l201ll and lAndrews et al.ll201ll . hereafter All). Better constraints on the 
disk structure can be obtained from reso lved images of the transitional disks. Using the Sub- 
millimeter Array (SMA) interferometer (JHo et al.ll2004j ). resolved imag es of transitional disks 
at sub-mm wavelength hav e provided direct detections of these cavities ( JAndrews fc Williams 
20071 : [Andrews et al.l 120091 ) . and measurements of their properties. Recently, All observed 
a sample of 12 nearby transitional disks (at a typical distance of ~140 pc). Combining both 
the SMA results and the SED, they fit detailed disk-|-cavity model for each object, and the 
cavity size (~15 — 70 AU) is determined with ~10% uncertainty. They concluded that large 
grains (up to ~mm-sized) inside the cavity are depleted by at least a factor of 10 to 100 
(the "depletion" in this work is relative to a "background" value extrapolated from the outer 
disk). Under the assumption that the surface density of the disk is described by their model, 
the infrared spectral fitting demands that the small grains (micron-sized and smaller) inside 
the cavity to be heavily depleted by a factor of ~10^. 

Recently, most objects in this sample have been observed by the Subaru High-Contrast 
Coronographic Imager for Adaptive Optics (HiCIAO) at NIR bands, as part of the th e 
Strategic Explorations of Exoplanets and Disks with Subaru project, SEEDS, (JTamurall2009l ). 
SEEDS is capable of producing polarized intensity (PI) images of disks, which greatly en- 
hances our ability to probe disk structure (especially at the inner part) by utilizing the fact 
that the central source is usually not polarized, so that the stellar residual in PI images 



^The objects with a small dust bel t left at the center have been prototyped as pre-transitional objects, 
such as LkCa 15 ( Espaillat et al.ll2008l ). but for simplicity we use the term transitional disk for both types. 
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is much smaller t han in full intensity (FI) images (jPerrin et al.l |2004J : iHinkley et al.l l2009l : 
Quanz et al.ll201ll ). 



The SEEDS results turned out to be a big surprise — in many cases the polarized 
NIR images do not show an inner cavity, despite the fact that the inner working angle 
of the images (the saturation radius or the coronagraph mask size, ipm ~ C'l ~ Gf!l5, or 
~15 — 25 AU at the distance to Taurus ~140 pc, see Section I2.3P is significantly smaller 
than the cavity sizes inferred from sub-mm observations. High contrast features such as 
surface brightness excesses or deficits exist in some systems, but they are localized and do 
not appear to be central cavities. Instead, the image is smooth on large scales, and the 
azimuthally averaged surface brightness radial profile (or the profile along the major axis) 
increases inward smoothly until ipiri, without any abrupt break or jump at the cavity edge (the 
slope may change with radius in some systems). Examples include ROX 44 (M. Kuzuhara 
et al. 2012, in prep.), SR 21 (K. Follet te et al. 2012, in prep.), GM Aur (J. Hashimoto et 
al. 2012a, in prep.), and SAO 206462 (JMuto et al.ll2012l ): see also the sample statistics (J. 
Hashimoto et al. 2012b, in prep.). Some objects such as UX Tau A also do not show a cavity 
(R. Tanii et al. 2012, in prep.), however the inner working angle of their SEEDS images is 
too close to the cavity size, so the status of the cavity is less certain. We note that LkCa 15 
also does not exhibit a clear cavity in its PI imagery (J. Wisniews ki et al. 2012, in prep .), 
but does exhibit evidence of the wall of a cavity in its FI imagery ( iThalmann et al.l |2010[ ) . 



This apparent inconsistency between observations at different wavelengths reveals some- 
thing fundamental in the transitional disk structure, as these datasets probe different com- 
ponents of protoplanetary disks. At short wavelengths (i.e. NIR) where the disk is optically 
thick, the flux is dominated by the small dust (micron-sized or so) at the surface of the 
disk (where the stellar photons get absorbed or scattered), and is sensitive to the shape of 
the surface; at long wavelength (i.e. sub-mm), disks are generally optically thin, so the flux 
essentially probes the disk sur face density in big grain s (mm-sized or so), due to their large 
opacity at these wavelengths (j Williams fc Ciezall201l[ ). 



Combining all the three pieces of the puzzle together (SED, sub-mm observation, and 
NIR imaging), we propose a disk model that explains the signatures in all three observations 
simultaneously: the key point is that the spatial distributions of small and big dust are 
decoupled inside the cavity. In this model, a well defined cavity (several tens of AU in radius) 
with a sharp edge exists only in spatial distribution of the big dust and reproduces the central 
void in the sub-mm images, while no discontinuity is found for the spatial distribution of the 
small dust at the cavity edge. Inside the cavity, the surface density of the small dust does 
not increase inwardly as steeply as it does in the outer disk; instead it is roughly constant 
or declines closer to the star (while maintaining an overall smooth profile). In this way. 
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the inner region (sub-AU to a few AU) is heavily depleted in small dust, so that the model 
reproduces the NIR flux deficit in the SED (but still enough small dust surface density to 
efficiently scatter near-IR radiation). Modeling results show that the scattered light images 
for this continuous spatial distribution of the small dust appear smooth as well, with surface 
brightness steadily increasing inwardly, as seen in many of the SEEDS observations. 

The structure of this paper is as follows. In Section [2] we introduce the method that we 
use for the radiative transfer modeling. In Section [3] we give the main results on the scattered 
light images: first a general interpretation of the big picture through a theoretical perspective, 
followed by the modeling results of various disk+cavity models. We investigate the sub-mm 
properties of these models in Section HI and explore the degeneracy in the disk parameter 
space on their model SED in Section |3 We summarize the direct constraints put by the 
three observations on this transitional disk sample in Section El as well as the implications 
of our disk models. Our generic solution, which qualitatively explains the signatures in all 
the three observations, is summarized in Section [71 



2. Radiative Transfer Modeling 

In this section, we introduce the model setup in our radiative transfer calculations, and 
the post processing of the raw NIR polarized scattered light images which we perform in 
order to mimic the observations. The purpose of this modeling exercise is to "translate" 
various physical disk models to their corresponding NIR polarized scattered light images, 
sub-mm emission images, and SED, for comparison with observations. 



2.1. Model setup 

We use a modifie d version of the Monte C arlo radiative transfer code developed by 



Whitnev et al.l ( l2003al lbl). iRobitaille et al. 



(120061). and B. W hitney et al. 2012, in prep.; for 



the disk structure, we use All and IWhitney et al.l (J2003bl ) for references. The NIR images 



(this section) and SED (Sectional) are produced from simulations with 4 x 10^ photon packets, 
and for the sub-mm images (SectionHD we use 5 x 10^ photon packets. By varying the random 
seeds in the Monte Carlo simulations, we find the noise levels in both the radial profile of the 
convolved images (Section l2.3p and the SED to be < 0.5% in the range of interest. In our 
models, we construct an axisymmetric disk (assumed to be at ~140 pc) 200 AU in radius 
on a 600 x 200 grid in spherical coordinates {R,6), where R is in the radial direction and 
6 is in the poloidal direction {6 = 0° is the disk mid-plane). We include accretion energy 
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in the disk using the Shakura & Sunyaev a disk prescription (jWhitney et aLll2003al ). Disk 
accretion under the accretion rate assumed in our models below (several x10~^Mq yr^^) 
does not have a significant effect on the SED or the images (for simplicity, accretion energy 
from the inn er gas disk is assurn ed to be emitted with the stellar spectrum, but see also the 



treatment in 



Akeson et al.l 120051 ) . We model the entire disk with two components: a thick 



disk with small grains (~/im-sized and smaller, more of less pristine), and a thin disk with 
large (grown and settled) grains (up to ~mm-sized). Figure [1] shows the schematic surface 
density profile for both dust population. 



The parametrized vertical density profiles for both dust populations are taken to be 
Gaussian (i.e. p(z) = poe~^ ^'^^ , isothermal in the vertical direction z\ with scale heights /ib 
and /is being simple power laws h oc R^ (we use subscripts "s" and "b" to indicate the small 
and big dust throughout the paper, while quantities without subscripts "s" and "b" are for 
both dust populations). Following All, to qualitatively account for the possibility of settling 
of big grains, we fix /ib = 0.2 x h^ in most cases to simplify the models, unless indicated 
otherwise. Radially the disk is divided into two regions: an outer full disk from a cavity edge 
i?cav to 200 AU, and an inner cavity from the dust sublimation radius -Rgub to i?cav (-Rsub is 
determined self -consistently as where t he temperature reaches the sublimation temperature 
Tlub ~1600 K, iDuUemond et al.ll200ll . usually around 0.1—0.2 AU). At places in the disk 
where a large surface area of material is directly exposed to starlight, a thin layer of ma terial 
is superheated, and the local disk "puffs" up vertically ( IDuUemond &: DominikI l2004a] ) . To 
study this effect at the inner rim (-Rsub) oi' at the cavity wall, we adopt a treatment similar 
to All. In some models below we manually raise the scale height h at -Rsub or -Rcav by a 
certain factor from its "original" value, and let the puffed up h fall back to the underlying 
power law profile of h within ~0.1 AU a s ^-(.^^l'^-^ ^u) We note that these puffed up walls 



are vertical, which may not be realistic (jlsella fc Nattal 120051 ). 



For the surface density profile in the outer disk, we assume 

-Re 



So(i?) 



''Cav^(R, 



-R)IRa 



R 



(1) 



where S^av is the surface density at the cavity edge (normalized by the total disk mass). 
Re is a characteristic scaling length, and the gas-to-dust ratio is fixed at 100. Following 
All, we take 85% of the dust mass to be in large grains aX R > -Rcav For the inner disk 
(i.e. R < -Rcav) three surface density profiles have been explored: 

Rr 



Si(-R) = ((5cavSc 



"'cav^(iJ, 



'R)/Rc 



R 



(rising Si(-R)), 



Si(-R) = ^cavScav (fiat Si(i?)), and 
Si(-R) = ((5cavScav)-B — (declining Si(-R)), 



(2) 
(3) 
(4) 



and their names are based on their behavior when moving inward inside the cavity. We note 
that Equation ([2]) and ([1]) together form a single S(i?) scahng relation for the entire disk 
(with different normalization for the inner and outer parts), as in All. 

We define the depletion factor of the total dust inside the cavity as 



Sp»(i?)^ 



m = ^;j^, (5) 



where S["^^(-R) is found by extrapolating So(-R) from the outer disk, i.e. evaluating Equa- 
tion (^ at R < i?cav (or Equation ([2]) with (5cav = !)• In addition, we define Ss{R) and 5b (-R) 
as the cavity depletion factors for the small and big dust respectively as 

'^^ 0.15 X Sf"ii(i?) ' "^^ 0.85 X Sf"ii(i?) ' ^^ 

where 0.15 and 0.85 are the mass fractions of the small and big dust in the outer disk. 
We note that unlike previous models such as All, our cavity depletion factors are radius 
dependent (a constant 6s{R) or 6h{R) means a uniform depletion at all radii inside the 
cavity). Specifically, we define the depletion factor right inside the cavity edge as 

^cav = S{R = -Rcav — c), (^cav.s = 4(-R = -Rcav — c), (^cav.b = '^b(-R = -Rcav — c)- (7) 

With the same ^cav, different models with different Ei(i?) profiles (Equations (I2l)-(j4j)) have 
similar Ei(/2) (and S{R)) in the outer part of the cavity, but very different Ei(i?) (and S{R)) 
at the innermost part. Lastly, the mass averaged cavity depletion factor (S) is defined as 

/„"'" Ef"(R)27cRdR 

SMA observations have placed strong constraints on the spatial distribution of the big 
dust, while the constraints on the small grains from the SED are less certain, especially 
beyond i? ~ 10 AU. Based on this, we adopt the spatial distribution of big grains in All 
(i.e. Equation ([1]), and no big grains inside the cavity), and focus on the effect of the distri- 
bution of small grains inside the cavity. Therefore, the sub-mm properties of our models are 
similar to those of the models in All (Section H]), since large grains dominate the sub-mm 
emission. We tested models with non-zero depletion for the big grains, and found that they 
make no significant difference as long as their surface density is below the SMA upper limit. 
From now on we drop the explicit radius dependence indicator (R) from various quantities 
in most cases for simplicity. 



-9- 



2.2. Dust properties 



For the sma ll grains we try two models: the standard interstel lar medium (ISM) g rains 
( iKim et al.lll994l . ~micron-sized and smaller), and the model that ICotera et al.l ( 120011 ) em- 
ployed to reproduce the HH 30 NIR scattered light images, which are somewhat larger than 
the ISM grains (maxim size ~20 /im). These grains contain silicate, graphite, and amor- 
phous carbon, and their properties are plotted in Figure |2l The two grain models are similar 
to each other, and both are similar to the small grains model which All used in the outer 
disk and the cavity grains which All used inside the cavity and on the cavity wall. We note 
that for detailed modeling which aims at fitting specific objects, the model for the small 
grains needs to be turned for each individual object. For example, the strength and shape 
of the silicate featu res i ndicate different co nditions for the small grains in the inner disk 
( lAdame et al.l l201ll and iFurlan et al.l l201ll who also pointed out that the silicate features 
in transitional disks typically show that the grains in the inner disk are dominated by small 
amorphous silicate grains similar to ISM grains). However, since we do not ain a at fitting 



speci fic objects, we avoid tuning the small dust properties and assume ISM grains (JKim et al. 
19941 ) for the models shown below, to keep our models generalized and simple. 



For the large grains we try three different models, namely Models 1, 2, and 3 from 
Wood et al.l (120021 ). The properties of these models ar e plotted in Figure HI They adopt a 



power-law size distribution (i.e. as in iKim et al.lll994j ) with an exponential cutoff at large 
size, and the maxim size is ~1 mm. These grains are made of amorphous carbon and 
astronomical silicates, with solar abundances of carbon and silicon. These models cover a 
large parameter space, however we find that they hardly make any difference in the scattered 
light image and the IRS SED, due to their small scale height and their absence ins i de the 
cavity. For this reason we fix our big grains as described by Model 2 in IWood et al.l (120021 ) 
(which is similar to the model of the big grains in All). We note that small grains have much 
larger opacity than big grains at NIR, and it is the other way around at sub-mm (Figure |2]). 



2.3. Post processing of the scattered light images 



To obtain realistic images which can be directly compared to SEEDS observations, the 
raw NIR images of the entire disk+star system from the radiative transfer simulations need 
to be convolved with the point spread function (PSF) of the instrument. SEEDS can ob- 
tain both the FI and the PI images for any object, either with or without a coronagraph 
mask. The observation could be conducted in several different observational modes, includ- 



ing an gular differential irn aging (ADI. iMarois et al.ll2006l ). polariz ation differential irn aging 
(PDI. iHinkley et al.ll2009l ). and spectral differential imaging (SDI. iMarois et al.ll2000l ). For 
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a description of the instrument see iTamural (120091 ) and ISuzuki et al.l ( 120101 ) 



In this work, we produce both the narrow band 880 /im images and H band NIR images. 
While at 880 /im we produce the full intensity images, for the NIR scattered light images we 
focus on the PI images (produced in the PDI mode, both with and without a coronagraph 
mask). This is because (1) PDI is the dominate mode for this sample in SEEDS, and (2) 
it is more difficult to interpret FI (ADI) images since its reduction process partially or 
completely subtracts azimuthally symmetric structure. Other authors had t o synthesize and 



reduc e model data in order to test for t he existence of features like cavities ( jThalmann et al. 



20101) or spatially extended emission (IThalmann et al.l 1201 ll ). For examples of PDI data 
reduction and analysis, see iHashimoto et al.l (1201 if ). When observing with a mask, ip^^ in 
the PI images is the mask size (typically 0'.'15 in radius), and when observing without a 
mask, 'j/'in is determined by the saturation radius, which typically is ~0'.'l. 

To produce an image corresponding to observations made without a mask, we convolve 
the raw PI image of the entire system with an observed unsaturated HiCIAO H-hand PSF. 
The re solution of the PSF is ~0'.'05 (~1.2A/D for an 8-m telescope) and the Strehl ratio is 
~40% ( ISuzuki et al.ll2010l ). The integrated flux within a circle of radius 0'.'25 is ~80% of the 



total flux (~90% for a circle of radius O'.'S). We then carve out a circle at the center with O'.'l 
in radius to mimic the effect of saturation. We call this product the convolved unmasked PI 
image. To produce an image corresponding to observations with a mask, we first convolve 
the part of the raw PI image which is not blocked by the mask with the above PSF. We 
then convolve the central source by an observed PI coronagraph stellar residual map (the 
PSF under the coronagraph), and add this stellar residual to the disk images (the fiux from 
the inner part of the disk which is blocked by the mask, ~20 AU at ~140 pc, is added to 
the star). Lastly, we carve out a circle O'.'IS in radius from the center from the combined 
image to indicate the mask. We call this product the convolved masked PI image. We note 
that the stellar residual is needed to fully reproduce the observations, but in our sample the 
surface brightness of the stellar residual is generally well below the surface brightness of the 
disk at the radius of interest, so it doesn't affect the properties of the images much. 

In this study, the disk is assumed to be face-on in order to minimize the effect of the phase 
function in the scattering, so that we can focus on the effect of the disk structure. This is a 
good approximation since most objects in this IRS/SMA/Subaru sample have inclinations 
around ~25° (i.e. minor to major axis ratio ~0.9. An observational bias towards face-on 
objects may exist, since they are better at revealing the cavity). Additional information 
about the scattering properties of the dust could be gained from analyzing the detailed 
azimuthal profile of the scattered light in each individual system, which we defer to the 
future studies. To calculate the azimuthally averaged surface brightness profiles, we bin the 
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convolved images into a series of annuli O'.'OS in width (the typical spatial resolution), and 
measure the mean flux within each annulus. 



3. The NIR polarized scattered light images 

With the tools described above, we investigate what kinds of disk structure could si- 
multaneously reproduce the gross properties of all three kinds of observations described in 
Section [1] In this section, we first investigate the properties of the scattered light images 
from a semi-analytical theoretical point of view (Section I3.ip . then we present the model 
results from the Monte Carlo simulations (Section 13. 2p . 



3.1. Theoretical considerations 

In a single NIR band, when the (inner) disk is optically thick (i.e. not heavily depleted 
of the small grains), the scattering of the starlight can be approximated as happening on 
a scattering surface Zg where the optical depth between the star and surface is unity (the 
single scattering approximation). This surface is determined by both the disk scale height 
(particularly /3 in the simple vertically-isothermal models), and the radial profile of the 
surface density of the grains. The surface br i ghtness of the sc attered light /scat(-R) scales 



with radius as (jJang-Condell fc Sasselovll2003l : llnoue et al.l 120081 ) 



/scat (R) ~ (P^^ sin 7 (9) 

where L^ is the stellar luminosity at this wavelength, C is a geometrical scattering factor, p is 
the polarization coefficient for PI {p = 1 for FI), and 7 is the grazing angle (the angle between 
the impinging stellar radiation and the tangent of the scattering surface). We note that both 
C and p depend on azimuthal angle, inclination of the disk, and the scattering properties of the 
specific dust population responsible for scattering at the particular wavelength. However, if 
the disk is relatively face-on and not too fiared, they are nearly position independent, because 
the scattering angle is nearly a constant throughout the disk and the dust properties of the 
specific dust population do not change much with radius. 

The grazing angle is determined by the curvature of the scattering surface. In axisym- 
metric disks, assuming the surface density and scale height of the small dust to be smooth 
functions of radius, the grazing angle is also smooth with radius, and two extreme conditions 
can be constrained as follows: 
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1. For disks whose scattering surface is defin ed by a constant poloidal a ngle 6 (such as a 
constant opening angle disk), 7 ~ -R*/-R (IChiang &: GoldreichI 119971 . where i?^ is the 
radius of the star), so the brightness of the scattered light scales with R as 



-XR) oc R- 



(10) 



2. For flared disks (but not too flared, h < R), the grazing angle can be well approximated 



as 



sm 7 ~ 7 



dR 



z^ 
R 



'.IV. 



In this case, iMutol ( 1201 ll ) ex plicitly calculated the pos ition of the scattering surface 



Zs at various radii (see also IChiang fc GoldreichI 119971 1 . and found Zs = rjhs where 
the coefficient r^ ~ a few and is nearly a constant. Combine d with the fact that 
K/R oc -R^~^ (0 ^ 1.25 — 1.3 as typical values in irradiated disks, IChiang &: Goldreich 



19971 : iHartmann et al.l Il998l ) , we have the intensity of the scattered light scales with 
radius as 

/3eat(i?) CX R^-\ (12) 



Although the above calculations are under two extreme conditions (for complete flat 
or flared disks), and they are based on certain assumptions and the observed images have 
been smeared out by the instrument PSF, the radial proflles (azimuthally averaged, or along 
the major axis in inclined systems) of SEEDS scattered light images for many objects in 
this sample lie between Equations f lTO|) and flT2|) in the radius range of interest. In order to 
guide the eye and ease the comparison between modeling results and observations, we use 
the scaling relation 

hc.t{R) (x R-^-' (13) 

to represent typical observational results, and plot it on top of the radiative transfer results, 
which will be presented in Section I3l2] (with arbitrary normalization). 

On the other hand, an abrupt jump in the surface density or scale height profile of the 
small dust in the disk produces a jump in 7 at the corresponding position. The effect of this 
jump will be explored in Section 13.21 



3.2. Modeling results 

First, we present the simulated H band PI images for a face-on transitional disk with 
a uniformly heavily depleted cavity with Si as Equation ([2]) (Figure |3]), and the associated 
surface brightness radial profiles (the thick solid curves in Figure H]). In each figure, the 
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three panels show the raw image, the convolved unmasked image, and the convolved masked 
image, respectively. This model is motivated by the disk+cavity models in All; we therefore 
use parameters typical of those models. Experiments show that the peculiarities in each 
individual All disk model hardly affect the qualitative properties of the images and their 
radial profiles, as long as the cavity is large enough (> 0'.'2) and the disk is relatively face-on. 

This disk harbors a giant cavity at its center with i?cav = 42 AU (~0'.'3 at 140 pc). The 
disk has the same inwardly rising S scaling both inside and outside the cavity as Equations 
([T]) and ([2]). Outside the cavity both dust populations exist, with the big/small ratio as 
0.85/0.15. Inside the cavity there is no big dust (5b = 0), and the small dust is uniformly 
heavily depleted to 6s = 5cav,s ~ 10~^ {6 ~ 10~^). The surface density profiles for both 
dust populations can be found in Figure [TJ The disk has total mass of O.OIM© (gas-to-dust 
mass ratio 100), h/R = 0.075 at 100 AU with /3 = 1.15, Re = 40 AU, and accretion rate 
M = 5 X 10~^Mo yr~^. The central source is a 3 Rq, 2Mq, 5750 K G3 pre-main sequence 
star. We puff up the inner rim and the cavity wall by 100% and 200%, respectively. 

The most prominent features of this model in both the unmasked and masked PI images 
are the bright ring at -Rcav, and the surface brightness deficit inside the ring (i.e. the cavity). 
Correspondingly, the surface brightness profile increases inwardly in the outer disk, peaks 
around -Rcav, and then decreases sharply. This is very different from many SEEDS results 
(such as the examples mentioned in Section [1]), in which both the bright ring and the inner 
deficit are absent, and the surface brightness radial profile keeps increasing smoothly all the 
way from the outer disk to the inner working angle, as illustrated by the scaling relation flT5|) 
in Figure m 

This striking difference between models and observations suggests that the small dust 
cannot have such a large depletion at the cavity edge. Figure H] shows the effect of uniformly 
filling the cavity with small dust on the radial profile for unmasked images (left) and masked 
images (right), leaving the other model parameters fixed. As Sg gradually increases from 
~10~^ to 1, the deviation in the general shape between models and observations decreases. 
The (5s = 1 model corresponds to no depletion for small grains inside the cavity (i.e. a full 
small dust disk); this model agrees much better with the scattered light observations, despite 
a bump around the cavity edge produced by its puffed up wall (Section [3.2.31) . although this 
model fails to reproduce the transitional-disk-like SED (Section |5]). 

We note that this inconsistency between uniformly heavily depleted cavity models and 
observations is intrinsic and probably cannot be solved simply by assigning a high polariza- 
tion to the dust inside the cavity. The polarization fraction (PI/FI) in the convolved disk 
images of our mo dels ranges from ~0.3 to ~0.5, comparable to observations (for example 



Muto et al.ll2012l ). Even if we artificially increase the polarization fraction inside the cavity 
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by a factor of 10, by scaling up the cavity surface brightness in the raw PI images (which 
results in a ratio PI/FI greater than unity) while maintaining the outer disk unchanged, 
the convolved PI images produced by these uniformly heavily depleted cavity models still 
have a prominent cavity at their centers. We also note that the contrast of the cavity (the 
flux deficit inside -Rcav) and the strength of its edge (the brightness of the ring at -Rcav) are 
partially reduced in the convolved images compared with the raw images. This is due both 
to the convolution of the disk image with the telescope PSF, which naturally smooths out 
any sharp features in the raw images, and to the superimposed seeing halo from the bright 
innermost disk (especially for the unmasked images). 

In the rest of Section [3] we focus on the small grains inside the cavity while keeping 
the big grains absent, and study the effect of the parameters Si, ^cav.s, and the puffing up 
of the inner rim and cavity wall on the scattered light images. The convolved images for 
three representative models are shown in Figure [5l and the radial profiles for all models are 
shown in Figure |6l Each model below is varied from one standard model, which is shown as 
the top panel in Figure [5] and represented by the thick solid curve in all panels in Figure |6l 
This fiducial model has a fiat S^ with no discontinuity at the cavity edge for the small dust 
(i.e. (5cav,s = 1), and no puffed up rim or wall, but otherwise identical parameters to the 
models above. 



3.2.1. The effect of T.; 

First we study the effect of three surface density profiles inside the cavity, namely rising 
(Equation ([2])), fiat (Equation ([2])), and declining E; (Equation (jl])). The surface densities 
for these models are illustrated in Figure [TJ Panel (a) in Figure E] shows the effect of varying 
the surface density on the image radial profiles. Except for shifting the entire curve up 
and down, different Sj produce qualitatively very similar images and radial profiles, and all 
contain the gross features in many SEEDS observations (illustrated by the scaling relation 
f lT5]) ). This is due to both the fact that smooth surface density and scale height profiles 
yield a smooth scattering surface, and the effect of the PSF. We note that there is some 
coronagraph edge effect at the inner working angle in the masked images, which is caused by 
that the part of the disk just outside (but not inside) the mask is convolved with the PSF. 
This results in a narrow ring of flux deflcit just outside the mask. In general, the flux is 



trustable beyond about one FWHM of the PSF from ^-.^ (~ O'.'OS + O'.'IS = 0'.'2) (iMuto et al. 



2OI2I . some instrumental effects in observations may also affect the image quality within one 



FWHM from the mask edge as well). 
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3.2.2. The effect of 5, 



cav,s 



We investigate the effect of different (5cav,s with ffat Si ([3]). Panel (b) in Figures |6] shows 
the effect on the radial profile of the convolved images. When deviating from the fiducial 
model with (5cav,s = 1 (i-e. a continuous small dust disk at -Rcav), a bump around the cavity 
edge and a surface brightness deficit inside the cavity gradually emerge. We quantify this 
effect by measuring the relative flux deficit at 0'.'2 (|-Rcav) as a function of (5cav,s, the small- 
dust discontinuity at the cavity edge (subpanel in each plot). For each model we calculate 
the ratio of the flux at 0'.'2, /0.2, to /o.e, the flux at 0'.'6 (= 2_Rcav), normalized by /0.2//0.6 
in the fiducial, undepleted model. For the model with a 50% discontinuity in Eg at -Rcav 
(<^cav,s = 0.5), /0.2//0.6 is ~15% lower than in the fiducial model in the convolved unmasked 
images, and ~20% lower in the convolved masked images. The latter is larger because the 
mask suppresses the halo of the innermost disk, thus the relative flux deficit in the masked 
image is closer to its intrinsic value, i.e. the deficit in the raw images. In this sense the masked 
images are better at constraining the discontinuity of the small dust than the unmasked 
images. The middle row in Figure shows the model images for (5cav,s = 0.3 (illustrated 
by the thin solid curve in Figure [1]) . The edge of the cavity is quite prominent in the raw 
image, while it is somewhat smeared out but still visible in the convolved images. For many 
objects in this IRS/SMA/Subaru cross sample, the SEEDS images are grossly consistent 
with a continuous small dust disk at the -Rcav, while in some cases a small discontinuity may 
be tolerated. 

For the purpose of comparison with the observations, we now discuss the detectability of 
a finite surface density discontinuity at the cavity edge, given the sensitivity and noise level 
of the SEEDS data (the numerical noise level in our simulations is well below the noise level 
in the observations, see Section l?!T]) . In typical SEEDS observations with an integration time 
of several hundred seconds, the intrinsic Poisson noise of the surface brightness radial profile 
due to finite photon counts is usually a few tenth of one percent at radius of inter ests, smaller 



than the error introduced by the instrument and the data reduction process. iMuto et al. 



(I2OI2I ) estimated the local noise level of the surface brightness to be ~10% at i? = 0'.'5 for 
the SEEDS SAO 206462 PI images, which should be an upper limit for the noise level in 
the azimuthally averaged surface brightness (more pixels) and in the inner region of the 
images (brighter). If this is the typical value for the instrument, then our modeling results 
indicate that SEEDS surface brightness measurements should be able to put relatively tight 
constraints on the surface density discontinuity for the small dust at -Rcav For example, 
SEEDS should be able to distinguish a disk of small dust continuous at -Rcav from a disk of 
small dust with a 50% density drop at -Rcav 

Thus, a lower limit on the small dust depletion factor at the cavity edge can be deduced 
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from detailed modeling for each individual object. This lower limit is likely to be higher 
than the upper limit from SMA on the depletion factor of the big dust (< 0.1 — 0.01), for 
objects with relatively smooth radial profiles. If this is confirmed for some objects in which 
the two limits are both well determined, it means that the density distribution of the small 
dust needs to somehow decouple from the big dust at the cavity edge. We will come back to 
this point in Section [61 



3.2.3. The effect of the puffed up inner rim and cavity wall 

Lastly, we explore the effects of the puffed up inner rim and the cavity wall on the 
images, by comparing the fiducial model (no puffing up anywhere) to a model with the inner 
rim puffed up by 100%, and another model with the cavity wall puffed up by 200%. Panel 
(c) in Figure [6] shows the effects. Puffing up the inner rim has little effect on the image, 
while the puffed up wall produces a bump at the cavity edge, similar to the effect of a gap 
edge. Images for the model with the puffed up wall are shown in Figure O (bottom row), 
where the wall is prominent in the raw image while been somewhat smeared out but still 
visible in the convolved images. Scattered light images should be able to constrain the wall 
for individual objects. Without digging deeply into this issue, we simply note here that the 
SEEDS images for many objects in this sample are consistent with no or only a small puffed 
up wall. 



4. The sub-mm properties and images 

4.1. The sub-mm intensity profile 

For all models shown in Section [3], the disk is optically thin in the vertical direction at 
880 /im (r ~ Sk;^ < 1, where ^i^ is the opacity per gram at z/). In this case, the intensity 
I^{R) at the surface of the disk at a given radius may be expressed as: 

/oo 
{Bu,TUz)i^u,hph{z) + B^^T4z)i<'u,sps{z)) dz (14) 

'OO 

where v = 341 GHz at 880 /zm, B^^Tiz) is the Planck function at T{z) (the temperature at 
2;), and p{z) is the vertical density distribution (T and p depend on R as well). Outside the 
cavity, the big dust dominates the 880 /im emission, since the big dust is much more efficient 
at emitting at ~880 /im than the small dust {Ku,h/i^u,s ^ 30 at these wavelengths). Inside the 
cavity there is no big dust by our assumption, so the sub-mm emission comes only from the 
small dust. The thick curves in Figure [7] show the 880 /im intensity as a function of radius 
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for the models in Section 13.2.11 (i.e. continuous small dust disks with rising, flat, or declining 
Si). In other words, this is an 880 /im version of the surface brightness radial profile for the 
"raw" image, without being processed by a synthesized beam dimension (the SMA version 
of the "PSF"). The bottom thin dashed curve is for the uniformly heavily depleted cavity 
model as in Figure |3] (4 = 5cav,s ~ 10~^, rising Ei) and no big dust inside the cavity, which 
represents the sub-mm behavior of the models for most systems in All. 

Since the d isk is roughly isothermal in the vertical direction near the mid-plane where 



most dust lies ( IChiang fc GoldreichI 119971 ) . T{z) may be approximated by the mid-plane 
temperature Tmid- The Planck function can be approximated as -Bv,t ~ 2v'^kT /(? at this 
wavelength due to hv <^ kT (880 /xm ~16 K, marginally true in the very outer part of 
the disk). In addition, since the two dust populations have roughly the same mid-plane 
temperature, but very different opacities {Ku,h ^ i^u,s), Equation flTl]) can be simplified to 
Iu{R) oc Tmid^i^^bSb for R > _Rcav and Iu{R) oc T^id^^^s^s for R < Rcav For our continuous 
small disk models with a complete cavity for the big dust, the intensity at 880 /xm drops by 
~2.5 orders of magnitude when moving from outside (/^(-Rcav + e)) to inside {Iu{Rca.v — e)) 
the cavity edge, due to both the higher opacity of the big dust and the fact that big dust 
dominates the mass at i? > i?cav Inside the cavity, the intensity (now exclusively from the 
small dust) is determined b y the factor TminS. For an i rradiated disk T^id increases inwardly. 



typically as Tmid oc R ^/^ (IChiang &: GoldreichI 119971 ). while Si in our models could have 



various radial dependencies (Equations (l2])-(jl])). In the fiat Si models (Equation ([3])), the 
intensity inside the cavity roughly scales with R as Ii,{R) oc i?^^/^ — the same as {Tmid(-R)} 
— and is 1.5 orders of magnitude lower than Ii,{Rca.v + e) in the innermost disk (around the 
sublimation radius). On the other hand, if there is no depletion of the small dust anywhere 
inside the cavity (i.e. the rising Si with 5s = 1 case), the intensity at the center (thick dashed 
curve) can exceed I^{Rcsv + e) by two orders of magnitude. 

In All, for 880 /im images of models with no big dust inside the cavity, the residual 
emission near the disk center roughly traces the quantity TSiK. In most cases. All found 
those residuals to be below the noise floor. Due to the sensitivity limit, the constraint on 
TSiK inside the cavity is relatively weak; nevertheless. All were able to put an upper limit 
equivalent to 5b ^ 0.01 — 0.1 for the mm-sized dust (with exceptions such as LkCa 15). Here 
we use a mock disk model to mimic this constraint. The top thin dashed curve in Figure [7| is 
from a model with uniform depletion factors 6s = S^ = 0.01 for both dust populations inside 
the cavity (so the entire disk has the same dust composition everywhere). The result shows 
that, qualitatively, various models with a continuous small dust disk and a complete cavity 
for the big dust are all formally below this mock SMA limit, though a quantitative fitting 
of the visibility curve is needed to constrain 6s and 6ca.v,s, in terms of upper limits, on an 
object by object basis. This may line up with another SED-based constraint on the amount 
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of small dust in the innermost disk, as we will discuss in the Section |5l 

4.2. The 880 fim images 

While the intensity discussion qualitatively demonstrates the sub-mm properties of the 
disks, Figure m shows the narrow band images at 880 fira for two disk models. The top row is 
from the model which produces Figure [3] (also the bottom thin dashed curve in Figure [7] and 
the left panel in Figured]), which is an All style model with a uniformly heavily depleted 
cavity with rising Si, 6s = 5cav,s ~ 10~^, and no big dust inside the cavity. The bottom row is 
from the fiducial model in Section \^72\ (which produces the top row in Figure O and the thick 
solid curve in Figure [7] and in the left panel in Figured]), which has a continuous distribution 
for the small dust with flat Si and Scs,v,s = 1, and no big dust at i? < i?cav as well. The 
panels are the raw images from the radiative transfer simulations (left), images convolved 
by a Gaussian profile with resolution ~0'.'3 (middle, to mimic the SMA observations. All) 
and ~0'.'l (right, to mimic future ALMA observations (Section 16. 3p . 

Both models reproduce the characteristic features in the SMA images of this transitional 
disk sample: a bright ring at the cavity edge and a flux deficit inside, agree with the semi- 
analytical analysis in Section UT] but very different NIK scattered light images. The intrinsic 
reason for this apparent inconsistency is, as we discussed above, that big and small dust 
dominate the sub-mm and NIR signals in our models, respectively. Thus two disks can have 
similar images at one of the two wavelengths but very different images at the other, if they 
share similar spatial distributions for one of the dust populations but not the other. 

Lastly, we comment on the effect of big to small dust ratio, which is fixed in this work 
as 0.85/0.15 to simplify the model (s ee the discussion of d epletion of the small dust in the 



surface layer of protoplanetary disks, iD'Alessio et al.ll2006l ). The scattering comes from the 



disk surface and is determined by the grazing angle, which only weakly depends on the 
small dust surface density, if it is continuous and smooth (Section [3.2. ip . Changing the mass 
fraction of the big dust in the outer disk from 0.85 to 0.95 in our (5cav,s = 1 models (effectively 
a factor of 3 drop in surface density of the small dust everywhere) introduces a ~20% drop 
in the surface brightness of the scattered light images, but a factor of 3 drop in the cavity 
880 /xm intensity {Iu{R) oc Ss at i? < -Rcav)- Lastly, we note that since the big-to-small dust 
sub-mm emission ratio is oc fi;;^^bSb/Kjy,sSs, the small grains must contain more than 90% of 
the total dust mass to dominate the sub-mm emission. 
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5. The transitional-disk-like SED 

In this section, we explore the parameter degeneracy in reproducing the transitional- 
disk-like SED with their distinctive NIR-MIR dips. SED fitting (particularly of the IRS 
spectrum) can only provide constraints on the spatial distribution of the small dust within 
a few or a few tens AU from the center, and it contains strong degeneracy in the parameter 
space (All). Below, we show that disk models with different cavity structures can produce 
roughly the same SED, containing the transitional disk signature, as long as their innermost 
parts are modestly depleted (by a factor of ~1000 or so). Except for the specifically men- 
tioned parameters, the other parameters of these models are the same as for the fiducial 
model in Section 13. 2( in particular there is no big dust inside the cavity. 



5.1. The degeneracy in SED fitting 

Figure [9] shows the SED for four disk models varied based on the fiducial model in 
Section 13.21 The model for the thick dashed curve has a uniformly heavily depleted cavity 
with rising S;, 4 ~ 10~^, and no big dust inside the cavity (illustrated by the thin dashed 
curve in the left panel of Figured]). The scale height profile has f3 = 1.15 and h/R = 0.085 
at 100 AU. The inner rim is puffed up by 100% and the outer wall is puffed up by 200%. The 
inclination is assumed to be 20°, Re = 15 AU, and i?cav = 36 AU. This model is motivated by 
the All disk+cavity structure. The full small dust disk model (the thin dashed curve) has 
otherwise identical properties but 5s = 1 (i.e. completely filled cavity for the small dust). 
The other two smooth small-dust disk models have much more massive inner disks with 
(^cav.s = 1 (i-e. a continuous small dust disk) and no puffed up inner rim or cavity wall. The 
solid curve model has fiat Si (Equation ([3]), illustrated by the thick solid curve in the left 
panel of Figured]), /3 = 1.33 and h/R = 0.08 at 100 AU. The dash-dotted curve model has 
declining Sj (Equation (jl]), illustrated by the dash-dotted curve in the left panel of Figured]), 
(3 = 1.25, and h/R = 0.078 at 100 AU. 

The two smooth small dust disk models with flat or declining Si produce qualitatively 
similar SED as the uniformly heavily depleted model (in particular, roughly diving to the 
same depth at NIR, and coming back to the same level at MIR, as the signature of transitional 
disks), despite the fact that they have very different structures inside the cavity. The minor 
differences in the strength of the silicate feature and the NIR flux could be reduced by tuning 
the small dust model and using a specifically designed scale height profile at the innermost 
part (around the sublimation radius or so). The main reasons for the similarity are: 

1. The depletion factor (or the surface density) at the innermost part (from Rsuh to 
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~1 AU or so). While the two smooth small dust disk models differ by ~5 orders of 
magnitude on the depletion factor (or the surface density) at the cavity edge from the 
uniformly heavily depleted model, the difference is much smaller at the innermost disk, 
where most of the NIR-MIR flux is produced. At the innermost disk, the small dust is 
depleted by ~3 orders of magnitude in the flat Si model, ~5 orders of magnitude for 
the declining Si model, and ~5 orders of magnitude in the uniformly heavily depleted 
model (with rising Si). On the other hand, the integrated depletion factor (4) for the 
small dust is ~0.3 for the fiat Si model, ~0.2 for the declining Si, and ~10^^ for the 
uniformly heavily depleted model, more in line with 6ca.v,s, because most of the mass is 
at the outer part of the cavity. We note that the total amount of small dust is not as 
important as its spatial distribution inside the cavity, and the amount of dust in the 
innermost part, in determining the NIR-MIR SED. 

2. The scale height of small grains hg at the innermost part. The two smooth small dust 
disk models are more flared than the uniformly heavily depleted model. While the 
three have roughly the same scale height outside the cavity, the difference increases 
inward. At 1 AU, hs for the uniformly heavily depleted model is 1.7x that of the fiat 
Si model and 2.4 x that of the declining Si model. 

3. The puffed up inner rim. The inner rim scale height is doubled in the heavily depleted 
model, which increases the NIR fiux and reduces the MIR fiux since the puffed up rim 
receives more stellar radiation and shadows the disk behind it. The puffed up inner 
rim is removed in the fiat or declining Si models. 

The surface density (or the depletion factor) and the scale height at the innermost 
part are considerably degenerate in producing the NIR to MIR fiux in the SED (All). In 
general, a disk which has a higher surface density and scale height at the innermost part 
and a puffed up inner rim intercepts more stellar radiation at small radii, and has more 
dust exposed at a high temperature, so it produces more NIR fiux. On the other hand, 
the shadowing effect cast by th e innermost disk on the outer disk causes less MIR emission 



( JDuUemond fc Dominikll2004al lbl). In this way, changes in some of these parameters could be 



largely compensated by the others so that the resulting SED are qualitatively similar. 

However, in order to reproduce the characteristic transitional disk SED, the value of the 
depletion factor inside the cavity cannot be too high. The increasing surface density at small 
radii would eventually wipe out the distinctive SED deficit, and the resulting SED evolves 
to a full-disk- like SED, as illustrated by the full small dust disk model in Figure M In our 
expe riments with not too flared (comparing with the canonical /3 ~ 1.25 — 1.3 in irradiated 



disk, IChiang fc GoldreichI llQQTt iHartmann et al.l Il998l ) , we flnd an upper limit on the order 
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of 10"^ for the depletion factor in the innermost part in our smooth disk models. We note 
that this limit depends on the detailed choices of the disk and cavity geometry, such as Re 
and -Rcav; and the big-to-small-dust ratio in the outer disk. 



5.2. Discussion of the disk model in producing the SED 

The inner rim is puffed up in All to intercept the starlight and to shadow material at 
larger radii. At a given radius, the scale height h of the gas disk scales as 

h ^ cjn, (15) 

where Q is the orbital frequency and Cg is the isothermal sound speed in the disk 

Cs ^ VkT/fi, (16) 

where T is the (mid-plane) temperature and /i is the average weight of the particles. If the 
scale heights of the dust and the gas are well coupled (e.g. for well-mixed-gas-dust models 
or a constant level of dust settling), tripling the rim scale height (not unusual in All) 
means an order of magnitude increase in its temperature. Due to the sudden change of the 
radial optical depth from ~0 to unity in a narrow transition region directly illuminated by 
the star, some puffin g up may be present, but probably not that significant. In addition. 



Isella &: Nattal ( 120051 ) pointed out that a reahstic puffed up rim has a curved edge (away 
from the star) instead of a straight vertical edge due to the dependence of Tgub on pressure, 
which further limits the ability of the rim to shadow the outer disk. Based on these reasons, 
we choose not to have the rim puffed up in our models, though we note that a relatively 
weak puffing up, as in the uniformly heavily depleted model here, does not make a major 
difference in the results. Similar idea applies to the puffed up wall as well. 

The typical /3 value assumed in All in this sample (/3 ~ 1. 15) is small compared with 



the canonical values f or irradiated disk models {(3 ~ 1.25 — 1.3, IChiang fc GoldreichI 119971 : 



Hartmann et al.lll998l ). This leads to that the temperature determined by the input scale 
height (Tinput, Equation (IT5l) and ( fT6l) ) may increase inwardly too steeply compared with 
the output mid-plane temperature calculated in the code (Toutput)- ^input at 100 AU in the 
three models are close to each other due to their similar scale height there, and all agree 
with Toutput (~30 K) within 20%. However, at 1 AU, while Tjnput in our smooth small dust 
disk models is close to Toutput (~220 K, within 30%), the input temperature in the uniformly 
heavily depleted model appears to be too high by a factor of ~3. 
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6. Discussion 

6.1. Direct constraints from observations and our generic model 

First, we review the direct, model- independent constraints on the disk structure which 
the three observations — the infrared SED, SMA sub-mm observations, and SEEDS NIR 
polarized scattered light imaging — put on many transitional disks in this cross sample: 

1. IRS reveals a distinctive dip in the spectra around 10 /xm, which indicates that the small 
dust (~/im-sized or so) in the inner part of the disk (from R^uh to several AU) must 
be moderately depleted. However, due to degeneracy in parameter space, the detailed 
inner disk structure is model dependent. Models with different cavity depletion factors. 
Si, and scale heights at the innermost disk could all reproduce the transitional disk 
signature. The IRS spectra are not very sensitive to the distribution of big dust. 

2. The SMA images show a sub-mm central cavity, which indicates that the big dust (mm- 
sized or so, responsible for the sub-mm emission) is heavily depleted inside the cavity. 
However, while the observations can effectively constrain the spatial distribution of the 
big dust outside the cavity, they can place only upper hmits on its total amount inside 
the cavity. SMA observations do not place strong constraints on the distribution of the 
small dust, though a weak upper limit for the amount of small dust inside the cavity 
may be determined based on the SMA noise level. 

3. SEEDS NIR polarized scattered light images are smooth on large scales, and have no 
clear signs of a central cavity. The radial profiles of many images increase inwardly 
all the way from the outer disk to the inner working angle without sudden jumps or 
changes of slope, indicating that the scattering surfaces and their shapes are smooth 
and continuous (outside ■j/'in). On the other hand, scattered light images are not very 
sensitive to the detailed surface density profiles and the total amount of small dust 
inside the cavity. The NIR images normally do not provide significant constraints on 
the distribution of the big dust. 

In this work, we propose a generic disk model which grossly explains all three observa- 
tions simultaneously. Previous models in the literature which assume a full outer disk and 
a uniformly heavily depleted inner cavity can reproduce (1) and (2), but fail at (3), because 
they also produce a cavity in the scattered light images, which contradicts the new SEEDS 
results. Through radiative transfer modeling, we find that qualitatively (3) is consistent 
with a smooth disk of small dust with little discontinuity in both surface density and scale 
height profile. Table [1] summarizes the key points in various models and compares their 
performances in these three observations. 
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Since we focus on generic disk models only which reproduce the gross features in obser- 
vations, and we do not try to match the details of specific objects, we more or less freeze many 
nonessential parameters in Sections [HE which do not qualitatively change the big picture 
for simplicity. The important ingredients include the dust properties (mostly for the small 
dust, both the size distribution and the composition). Eg, the big-to-small-dust ratio, and hi 
(both the absolute scale and (3). NIR scattered light images are able to provide constraints 
on some of these parameters (particularly h^ and /3), due to the dependence of the position 
and the shape of the disk surface on them. These parameters were not well constrained 
previously using sub-mm observations and SED due to strong degeneracies (All). 

We note that alternative models for explaining the scattered light images exist, but 
generally they require additional complications. As one example, if the small dust is not 
depleted in the outer part of the sub-mm cavity, but is heavily depleted inside a radius 
smaller than ip^^ (< 15 — 20 AU), then it is possible to fit all the three observations, in which 
case the small dust cavity does not reveal itself in the scattered light images due to its small 
size. However, in this case one needs to explain why different dust populations have different 
cavity sizes. Future scattered light imaging with even smaller ipi^^ may test this hypothesis. 

As another example, while we achieve a smooth scattering surface by having continuous 
surface density and scale height profiles for the small dust, it is possible to have the same 
result with discontinuities in both, but with just the right amount such that the combination 
of the two yields a scattering surface inside the cavity smoothly joining the outer disk. 
This may work if the cavity is optically thick (i.e. not heavily depleted), so that a well- 
defined scattering surface inside the cavity exists. Experiments show that for the models in 
Section 13. 2[ uniformly depleting the small dust by a factor of ~1000 and tripling the scale 
height inside the cavity would roughly make a smooth scattering surface. However, fine 
tuning is needed to eliminate the visible edge from a small mismatch in the two profiles. 
Also, the thicker cavity shadows the outer disk, and makes it much dimmer in scattered 
light (by about one order of magnitude). Lastly, without tuning on the scale height and/or 
surface density in the innermost part, this model produces too much NIR-MIR flux and too 
small flux at longer wavelengths in its SED, due to its big scale height at small R and the 
subsequent shadowing effect. 



6.2. The structure of the cavity in transitional disks 

There are several important conclusions that can be drawn based on our modeling of 
the transitional disks at different wavelengths. 



-24- 



First, as we discussed in Section |3. 2. 2[ for some objects the lower limit for the depletion 
of the small dust at the cavity edge (as constrained by the scattered hght images) is hkely 
to be above the upper hmit for the big dust constrained by the SMA, based on the modehng 
resuhs and the noise level in the two instruments. This essentially means that the small 
dust has to spatially decouple from the big dust at the cavity edge. This is the first time 
that this phenomenon has been associated with a uniform sample in a systematic manner. 
Detailed modeling of both images for individual objects is needed, particularly in order to 
determine how sharp the big dust cavity edge is from the sub-mm observations, to pin down 
the two limits and check if they are really not overlapping. While we defer this to future 
work, we note that having the big dust the same surface density as the small dust inside the 
cavity probably cannot reproduce the sub-mm images. Experiments show that even with 
our declining Sj (jl]), a fixed big/small dust ratio and a continuous surface density for both 
throughout the disk (i.e. Sca.v,s = (^cav.b = 1) produce a sub-mm cavity with a substantially 
extended edge, and the central flux deficit disappears in the smeared out image. If this 
is confirmed, it further leads to two possibilities: (a) whatever mechanism responsible for 
clearing the cavity have different efficiency for the small and big dust, or (b) there are other 
additional mechanisms which differentiate the small and big dust after the cavity clearing 
process. At the moment it is not clear which one of the two possibilities is more likely, and 
both need more thorough investigations. 

Second, as we argued in Section [5l in order to reproduce the distinctive NIR deficit in 
the transitional disk SED, an effective "upper limit" of 6s at the innermost region is required, 
which in experiments with our disk parameters is on the order of 10~^. This is far from the 
lower limit of 4 at the cavity edge (close to 1), constrained by the scattered light images. 
Together, the two limits indicate that the spatial distribution of the small grains is very 
different inside and outside the cavity — specifically. Eg tends to be fiat or even decrease 
inwardly inside the cavity. 

In addition, this implies that the gas-to-dust ratio needs to increase inwardly, given that 
most of these objects have non-trivial accretion rates (M ~ 10~® — lO~^M0/yr, All). For 
a steady Shakura & Sunyaev disk, the accretion rate M is related to the gas surface density 

where a is the Shakura-Sunyaev viscosity parameter, Q is the angular velocity of the disk 
rotation, and Cg is given by Equation flT6|) . At ~0.1 AU, equation f lT7|) predicts Sgas ~10^ g 
cm~^, assuming a temperature T ~ 10^ K, M ~ lO~^M0/yr, a ~ 0.01, and M ~ Mq as 
typical T Tauri values. This is very different from our upper limit of Sgas ~1 g cm~^ in 
the innermost disk, obtained assuming a fixed gas-to-dust ratio of 100 (the fiat Si models 
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in Figured]). To simultaneously have large Sgas but small Sdust in the innermost disk, the 
gas-to-dust ratio needs to increase substantially fr om the nominal value of 100 (by a factor 
of ~ 10^ in our models, echos with lZhu et al.ll201ll ). This could put constraints on the cavity 
depletion mechanism or dust growth and settling theory. 

At the moment, the mechanism(s) which are responsible for clearing these giant cav- 
ities inside transitional disks are not clear (see summary of the current situation in All). 
Regarding the applications of our model on this subject, we note two points here. The 
flat /declining surface density of the small grains inside the cavity in our models is consistent 
with the grain growth and settling argument, that the small grains in the inner disk are con- 
sumed at a faster rate due to higher growth ra te there (JDuUemond &: Dominikll2005l : iGaraud 
20071 : iBrauer et al.ll2008l : iBirnstiel et al.ll2010l ). Also, the so-called dust filtration mechanism 
seems promising for explaining why small dust b ut not big dust is present inside the cavity 
( jPaardekooper fc Mellemall2006t iRice et al.ll2006l ). since it could effectively trap the big dust 
at a pressure maximum in the disk but filter th rough the s mall dust. Particularly, com- 
bining the two (dust growth and dust filtration), IZhu et al.l (120121 ) proposed a transitional 
disk formation model from a theoretical point of view to explain the observations, and their 
predicted spatial distribution of both dust populations in the entire disk is well consistent 
with the ones here. 



6.3. Future observations 



Imaging is a very powerful tool for constraining the structure of protoplanetary disks 
and the spatial distribution of both the small and big dust, and there are many ongoing 
efforts aiming at improving our ability to resolve the disks. In the direction of optical-NIR 
imaging, updating existing coronagraph and Adaptive Optics (AO) systems, such as the new 



Corqnagraphic Extreme Adaptive Optics (SCExAO) system on Subaru (JMurakami et al. 
2010l : iGuyon et al.l |20 111 which could raise the Strehl ratio to ~0.9), are expected to achieve 
better performance and smaller tp[^ in the near future. 

To demonstrates the power of the optimal performances of these next generation instru- 
ments in the NIR imaging. Figure [10] shows the surface brightness radial profile of several 
masked i7-band disk images convolved from the same raw image by different PSF. Except 
the dotted curve, all the other curves are from the model corresponding to the middle row 
in Figure [S] which has a 0'.'3 radius cavity, fiat Si, and 5cav,s = 0.3 (a 70% drop in Eg at 
-Rcav)- We use three PSFs: the current HiCIAO PSF in H band, which could be roughly ap- 
proximated by a diffraction limited core of an 8-m telescope (resolution ~0'.'05) with a Strehl 
ratio of ~0.4 plus an extended halo; mock PSF I (to mimic SCExAO), which is composed 
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of a diffraction limited core of an 8-ni telescope with a Strelil ratio of 0.9, and an extended 
halo similar in shape (but fainter) as the current HiCIAO PSF; model PSF II (to mimic the 
next generation thirty-to- forty meter class telescopes), which is composed of a diffraction 
limited core of a 30-m telescope (resolution ~0'.'013) with a Strehl ratio of 0.7, and a similar 
halo as the previous two. Compared with the full small dust disk case, all the convolved 
images of the 5cav,s = 0.3 model shows a bump at i?cav and a relative flux deficit at R < Raw- 
However, from the current HiCIAO PSF to model PSF I and II, the contrast level of the 
cavity becomes higher and higher, and closer and closer to the raw image (which essentially 
has an infinite spatial resolution). With these next generation instruments, the transition of 
the spatial distribution of the small dust at the cavity edge will be better revealed. 

On the other hand, in radio astronomy, the Atacama Large Millimeter Array (ALMAcI) 
is expected to revolutionize the field, with its much better sensitivity level and exceptional 
spatial resolution (~0'.'l or better). As examples, the right panels in Figure [8] show images 
convolved by a Gaussian profile with resolution ~0'.'l, which mimic the ability of ALMA and 
show two prominent improvements over the images under the current SMA resolution (~0'.'3, 
middle panels). First, the edge of the cavity is much sharper in the mock ALMA images. 
This will make the constraint on the transition of the big dust distribution at the cavity edge 
much better. Second, while the weak emission signal in the bottom model is overwhelmed 
by the halo of the outer disk in the mock SMA image, resulting in that the bottom model 
is nearly indistinguishable from the top model (which essentially produces zero cavity sub- 
mm emission), the mock ALMA image successfully resolves the signal as an independent 
component from the outer disk, and separates the two models. This weak emission signal 
traces the spatial distribution of the dust (both populations) inside the cavity, which is the 
key in understanding the transitional disk structure. 

At this stage, the total number of objects which have been observed by all three survey- 
scale projects (using IRS/SMA/Subaru) is still small. Increasing the number in this multi- 
instrument cross sample will help clear the picture. In addition, future observations which 
produce high spatial resolution imag es at other wavelen gths, such as UV, optical, or other 
NIR bands (for example using HST, iGrady et al.ll2009l . or the future JWST), or using in- 
terferometer (such as the Astrometric and Phase- Referencing Astronomy project on Keck, 



Adelman et a 



20071: IPott et al.l l2009l. and AMBER system on Very Large Telescope In- 



Petrov et al.l 120071 : iTatuUi et al.l 120071 ) should also be able to provide useful 



terferometer, 

constraints on the disk properties. 



^http://www. alniaobservatory.org/ 
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7. Summary 

We summarize this paper by coming back to the question which we raised at the begin- 
ning: what kind of disk structure is consistent with and is able to reproduce the characteristic 
signatures in all three observations of transitional protoplanetary disks: a high contrast cav- 
ity in sub-mm images by SMA, a NIR deficit in SED by Spitzer IRS, and a smooth radial 
profile in NIR polarized scattered light images by Subaru HiCIAO. We propose one generic 
solution for this problem, which is feasible but by no means unique. The key points are: 

1. A cavity with a sharp edge in the density distribution of big grains (up to ~mm-sized) 
and with a depletion factor of at least 0.1-0.01 inside is needed to reproduce the SMA 
sub-mm images, as pointed out by All. 

2. Right inside the cavity edge (~15-70 AU), the surface density for the small dust 
(~niicron-sized and smaller) does not have a big sudden (downward) jump (A small 
discontinuity may exist). The SEEDS NIR scattered light images, which typically 
detect the disk at R >15 AU (the inner working angle in SEEDS), generally require 
continuous/smooth profiles for the surface density and scale height of the small dust. 

3. The small dust in the innermost region (i.e. within a few AU, on a scale smaller 
than measured by SEEDS) has to be moderately depleted in order to produce the 
transitional-disk-like SED, assuming the disk is not too flared, but the exact depletion 
factor is uncertain and model dependent. 



As we discussed in Section 16. 2[ combining all the above points, our model suggests that 
the spatial distributions of the big and small dust are decoupled mside the cavity (particularly 
at the cavity edge). Also, our model argues that the surface density of the small dust inside 
the cavity is flat or decreases with radius, consistent with the predictions in dust growth 
models. Combined with the accretion rate measurement of these objects, it further implies 
that the gas-to-dust ratio increases inwardly inside the cavity of transitional disks. 
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Fig. 1. — Various surface density profiles of the small dust (up to ~/im-sized, left) and big 
dust (up to ~mm-sized, right) in our model (gas-to-dust ratio is assumed to be 100:1 in this 
work). The big grains contain 85% of the dust mass outside the cavity, and are completely 
absent inside. The thick curves at -R < -Rcav in the left panel indicate three different Si (rising, 
fiat, and declining, as Equation (l2])-(j4j), with names indicating their behavior when moving 
inward inside the cavity) for the small dust, all continous at the cavity edge (5cav,s = !)• The 
bottom thin dashed curve illustrates the surface density profile in typical All models, which 
have a uniformly heavily depleted cavity with 4 ~ 10~^, and the thin solid curve represents 
a slight small-grain depletion at the cavity edge (5cav,s = 0.3) with a fiat Ei. The vertical 
dotted curves indicate the cavity edge and typical SEEDS inner working angle (~0.1 — 0.15" 
or ~15 — 20 AU at Taurus). The labels on the axes are for Section |3] and HI 
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Table 1: Model Comparison 





Our model 


Uniformly heavily depleted cavity 


Full small dust 


Si,s 


Equation ©-(H 


Equation (^ 


Equation (^j 


Os at -fXcav (Ocav.s) 


1 


-10-5 


1 


6s at -Rsub 


-10-5 - 10-3 


-10-5 


1 


Si^b (or 5b) 











hs at -Rsub 


thinner 


thicker 


— 


puffed up inner rim 


no or little 


little to heavy 


— 


observation*^ : SED 


Y 


Y 


N 


observation: sub-mm 


Y 


Y 


Y 


observation: NIR 


Y 


N 


Y 



"Whether or not these models can reproduce the signatures in various observations: the infrared deficit around 
10 microns in SED, the central cavity in sub-mm images, and the smooth radial profile of NIR images. 
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Fig. 2. — Dust properties of the grain models we use, showing opacity k, (a), albedo u (b), 
average cosine scattering angle g (c), and maximum polarization Pi (d) (which typically 
occurs at a 90° scattering angle). The opacity is the dust+gas opacity assuming a gas-to- 
dust ratio of 100. Thick curves are for the small dust models, which dominate the scattered 
light due to their large scale height and NIR opacity. Thin curves are for the big dust models, 
which dominate the sub-mm emission due to their high opacity at these wavelengths. The 
solid curves are the fiducial choices in this work. 
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Fig. 3. — Simulated H band scattered light PI images for a face-on disk with a uniformly 
heavily depleted (4 ~ 10^^) cavity (see Section [3l2] for details), displayed with logarithmic 
stretches. The left panel shows the raw disk image from Whitney's radiative transfer simula- 
tion, the middle panel shows the unmasked image convolved with the HiCIAO PSF, and the 
right panel shows the convolved masked image (see Section 12.31 for details on the HiCIAO 
PSF). The cavity is 42 AU in radius (~0'.'3 at 140 pc), and the inner working angles are 
~0'.'l without a coronagraph, and ~0'.'15 with a coronagraph (the hole in the middle and 
right panels, respectively). The central surface brightness deficit and the bright ring at the 
cavity edge are prominent under both observation modes, which is qualitatively inconsistent 
with many SEEDS transitional disks observations. This indicates that the cavity is not 
heavily depleted in small dust. 
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Fig. 4. — Surface brightness radial profiles of the raw images (a), convolved unmasked PI 
images (b), and convolved masked PI images (c), assuming a cavity completely devoid of 
large dust and uniformly depleted in small dust. The model with 5^ = 10~^ (thick solid 
curves) corresponds to Figure [3] and typical models in All, while the (5s = 1 model (dotted 
curves) has no depletion in small grains at the cavity edge. The thin solid curves labeled as 
Equation flT3l) represent the scaling relation typical in SEEDS (with arbitrary normalization). 
The basic features of the uniformly heavily depleted models — an increase in intensity at 
the cavity edge and a drop inside — are inconsistent with typical SEEDS results, which 
smoothly increase in intensity up to the inner working angle. 
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Fig. 5. — Simulated H band scattered light PI images for three face-on disks with a 42 AU 
cavity (~0'.'3 at 140 pc), which is devoid of big dust and filled with a uniform surface density 
of small dust (flat Si, Equation ([3])), displayed with logarithmic stretches. Panels (from left 
to right) are raw images, convolved unmasked images, and convolved masked images (the 
same as in Figure |3]). The top and middle rows show cavities undepleted and 70% depleted 
in small dust at the edge ((5cav,s = 1 and 5cav,s = 0.3, respectively), corresponding to the solid 
curves in Figure |3l The bottom row shows a cavity undepleted in small dust at the edge 
but with the cavity wall puffed up by 200%. For detailed model parameters see Section [3^ 
Images in the top row are qualitatively similar to many SEEDS transitional disk observations 
in this sample, which are smooth on large scales and lack abrupt breaks or bumps at the 
cavity edge. Also, results show that both a modest surface density discontinuity (a 70% 
drop, middle row) and a puffed up cavity wall (bottom row) could well be visible. 
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Fig. 6. — Surface brightness radial profiles of the convolved unmasked (left) and masked 
(right) PI images for the models in Section |3T2| showing the effect of Si (a), 5cav,s (b), and 
the puffed up inner rim and cavity wall (c). The subpanel in panel (b) shows the relative 
flux deflcit at 0'.'2 (|-Rcav) as a function of 5cav,s (Section |3. 2. 2p . All models have no big dust 
inside the cavity. The curves labled "flducial" are the same in all panels, which is from a 
model with flat Ej, a continuous small dust surface density proflle ((5cav,s = 1), and no puffed 
up cavity wall (the top row in Figure [5]). Results show that both a modest depletion in the 
small grains at -Rcav and a puffed up cavity wall can add a bump to the originally smooth 
radial proflle, while various forms of Si with no discontinuity at i?cav produce similar surface 
brightness radial proflle. Typical SEEDS results, as summarized in Section [1] and illustrated 
by the scaling relation f lT^ (with arbitrary normalization), is consistent with smooth small 
dust disks continuous in both surface density and scale height. The apparent inward surface 
brightness decrease near ipij^ in the masked images is a coronagraph edge effect (Section l3.2.ip . 
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Fig. 7. — The 880 /iin intensity Iy{R) profile (Equation flT4|) ) for disk models witli various 
inner disk structures. Tlie model for the top thin dashed curve has a uniformly depleted 
cavity with 5^ = 5^ = 0.01 (i.e. the same dust composition as the outer disk), in order to 
mimic a "mock" SMA upper limit on dust inside the cavity (All), while all the other models 
have no big dust (~mm-sized) inside the cavity. The three thick curves are from models 
with different Ej and no discontinuity at the cavity edge, corresponding to the same curves 
in Figure [1] and panel (a) in Figure [61 The bottom thin dashed curve is from a model with 
a uniformly heavily depleted cavity {5s ~ 10~^). See Section WA\ for details of the models. 
Qualitatively, the disk models with no big dust inside the cavity all produce a deficit in the 
intensity inside the cavity, since the small dust is very inefficient at sub-mm emission, and 
all are consistent with the mock upper limit set by the SMA results (All). 
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Fig. 8. — Narrow band 880/im images for two disk models with no big dust inside the 
cavity, displayed with logarithmic stretches. Panels are the raw images from the radiative 
transfer simulations (left), the processed imaged convolved with a Gaussian profile with 
resolution 0'.'3 (middle), to mimic SMA observations (All), and with resolution O'.'l (right), 
to mimic ALMA observations. The top images are from the uniformly heavily depleted 
model in Figure [3] (rising Si with 6s ~ 10~^, also the thin dashed curve inland left panel 
in Figured]). The bottom ones are from the same model as the top row in Figure [5] (also the 
thick solid curve in Figure [7] and left panel in Figure [1]), which has a continuous distribution 
for the small dust ((5cav,s = 1) with flat Sj (Equation ([3])). See Section lO] for details of the 
models. The two models produce similar sub-mm images but very different NIR scattered 
light images (Figure [3] and [5l). The bottom right panel shows that with its exceptional spatial 
resolution (~0'.'l or better), ALMA may reveal a sharper edge of the cavity and resolve the 
small sub-mm emission inside the cavity, which are the keys in a better understanding of the 
spatial distribution of the grains. 
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Fig. 9. — The SED for four disk+cavity models. The uniformly heavily depleted model is 
motivated by All, which has an inner cavity with 6s ~ 10~^, rising S^, and puffed up inner 
rim and cavity wall. The full small dust disk model (the thin dashed curve) has otherwise 
identical properties but with the small dust cavity being uniformly filled up (4 = 1). The 
other two smooth small dust disk models have (5cav,s = 1 (no discontinuity at the cavity edge) 
and no puffed up inner rim or wall. See Section [H?T] for detailed parameters in these models. 
The three thick curves show the parameter degeneracy in producing the SED — diving to 
roughly the same depth at NIR and coming back to the same level at MIR, as the signature 
of transitional disks. The full small dust disk model shows that increasing surface density 
at small radii would eventually wipe out the distinctive deficit in the NIR, and the resulting 
SED gradually evolves to a full-disk-like SED. 
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Fig. 10. — Examples of the power of the optimal performances of the next generation instru- 
ments in NIR imaging (such as the SCExAO on Subaru and thirty meter class telescopes). 
The thick solid, dashed, and dot-dashed curves show the surface brightness radial profile of 
several H band masked disk images convolved from the same raw image by several different 
PSF (indicated by the legend, see Section [673] for details). The base model corresponds to 
the middle row in Figure |5l which has a 0'.'3 radius cavity, flat Sj, and 5cav,s = 0.3 (a 70% 
drop in Eg at -Rcav)- No coronagraph stellar residual (see Section 1231) is added in order to 
isolate the effect of the PSF, and we bin the images into annulus 0'.'02 in width for better 
illustration. For comparision, the thin solid line is for the raw (unconvolved) image, and the 
dotted line is from the corresponding full small dust disk model (i.e. the top row in Figure |5]). 
These examples show that with the ability of next generation instruments expected in the 
next decade or so, the transition of the spatial distribution of the small dust at the cavity 
edge could be constrained much better. 



